//static void printConfigBackup(caja *bx, parametros param)
//{
//	char name[50];
//	sprintf(name, "%s","system_bk.dat");
//	ofstream MiArchivo (name,ios::out);
//	MiArchivo.precision(4);
//	if (MiArchivo.is_open()) {
//	    MiArchivo.setf(ios::fixed);
//		MiArchivo<<"Parametros de la simulacion"<<endl;
//		MiArchivo.fill('*');
//		MiArchivo << setw(122) << ""<<endl;
//		MiArchivo.fill(' ');
//		MiArchivo<<"Npart"<<setw(10)<<"phi"<<setw(10)<<"rcut"<<setw(10)<<"Dt"<<setw(10)<<"restart"<<setw(10)<<"ncorr"
//			<<setw(10)<<"trelaj1"<<setw(10)<<"bx"<<setw(10)<<"by"<<setw(10)<<"bz"<<setw(10)<<"tiempo"<<setw(10)<<"paso"<<endl;
//		MiArchivo.fill('-');
//		MiArchivo << setw(122) << ""<<endl;
//		MiArchivo.fill(' ');
//		MiArchivo<<bx->npart<<setw(10)<<bx->phi<<setw(10)<<param.rcut<<setw(10)<<param.Dt<<setw(10)<<param.rest<<setw(10)<<param.ncorr
//			<<setw(10)<<param.trelaj1<<setw(10)<<bx->lbox[0]<<setw(10)<<bx->lbox[1]<<setw(10)<<bx->lbox[2]<<setw(10)<<param.Dt*bx->nacc<<setw(10)<<bx->nacc<<setw(10)<<endl;
//	}
//
//	MiArchivo<<endl;
//	MiArchivo<<"Posiciones y velocidades"<<endl;
//	MiArchivo.fill('*');
//	MiArchivo << setw(122) << ""<<endl;
//	MiArchivo.fill(' ');
//	MiArchivo<<"rx"<<setw(20)<<"ry"<<setw(20)<<"rz"<<setw(20)<<"vx"<<setw(20)<<"vy"<<setw(20)<<"vz"<<endl;
//	MiArchivo.fill('-');
//	//MiArchivo << setw(122) << ""<<endl;
//	MiArchivo.fill(' ');
//	for(int i=0;i<bx->npart;i++) {
//	    MiArchivo<<bx->part[i].pos[0]<<setw(20)<<bx->part[i].pos[1]<<setw(20)<<bx->part[i].pos[2]<<setw(20)<<bx->part[i].vel[0]<<setw(20)<<bx->part[i].vel[1]<<setw(20)<<bx->part[i].vel[2];
////	    #ifdef POLYDISPERSE
////            MiArchivo<<setw(20)<<(bx)->part[i].radio;
//        //MiArchivo<<1<<endl;
////	    #endif
//	    //MiArchivo.fill(' ');
//        MiArchivo<<endl;
//	}
//
//	MiArchivo.close();
//}

static void printConfigBackup(caja *bx, parametros param)
{
	char name[50];
	sprintf(name, "%s","system_bk.dat");
	ofstream MiArchivo (name,ios::out);
	if (MiArchivo.is_open()) {
		MiArchivo<<bx->npart<<"  "<<bx->phi<<"  "<<param.rcut<<"  "<<param.Dt<<"  "<<param.rest<<"  "<<param.ncorr
			<<"  "<<param.trelaj1<<"  "<<bx->lbox[0]<<"  "<<bx->lbox[1]<<"  "<<bx->lbox[2]<<"  "<<param.Dt*bx->nacc
			<<"  "<<bx->nacc<<"  "<<endl;
	}

	for(int i=0;i<bx->npart;i++) {
	    MiArchivo<<bx->part[i].pos[0]<<"  "<<bx->part[i].pos[1]<<"  "<<bx->part[i].pos[2]<<"  "
	    <<bx->part[i].vel[0]<<"  "<<bx->part[i].vel[1]<<"  "<<bx->part[i].vel[2];
	    #ifdef POLYDISPERSE
            MiArchivo<<"  "<<(bx)->part[i].radio;
        //MiArchivo<<1<<endl;
	    #endif
	    //MiArchivo.fill(' ');
        MiArchivo<<endl;
	}

	MiArchivo.close();
}


static void printShearStressTensor(caja *box, parametros param)
{
	char name[50];
	sprintf(name, "%s", "stressTensor.dat");
	ofstream MiArchivo2 (name,ios::out|ios::app);
	if (MiArchivo2.is_open())
	{
		MiArchivo2 << box->nacc*param.Dt<<"  "<<box->svxy+box->spxy<<"  "<<box->svxz+box->spxz
		<<"  "<<box->svyz+box->spyz<<"  "<<box->svzz+box->spzz<<endl;
        MiArchivo2.close();
	}
}

static void printXYZfile(caja *box)
{
	char name[50];
	sprintf(name, "%s", "representacion.dat");
	ofstream MiArchivo (name,ios::out|ios::app);
	MiArchivo.precision(5);
	if (MiArchivo.is_open())
	{
		MiArchivo << box->npart<<endl;
		MiArchivo << box->lbox[0]<<"  "<<box->lbox[1]<<"  "<<box->lbox[2]<<"  "<<box->nacc<<endl;
		for(int i=0;i<box->npart;i++)
		{
			MiArchivo<<box->part[i].pos[0]<<"  "<<setw(8)<<box->part[i].pos[1]<<"  "<<setw(8)<<box->part[i].pos[2]<<"  "<<setw(8)<<box->part[i].radio<<endl;
		}
		MiArchivo.close();
	}
}
